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We identify and quantify systematic effects not accounted for in two previous measurements of tiie 
a-a relative-energy distribution in the /3 decay of *B, wiiidi can explain the apparent disagreement 
with respect to two newer measurements. This settles a current dispute concerning the shape of the 
*B neutrino spectrum of importance to solar-neutrino studies. 
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Introduction. New solar-neutrino observatories arc 
allowing us to test our understanding of the Sun and 
of neutrinos themselves with increasing precision. Neu- 
trinos emitted in the (3 decay of play a central role 
in these endeavours because they dominate the solar- 
neutrino spectrum above approximately 2 MeV. Precise 
and accurate knowledge of the laboratory '^B neutrino 
spectrum is necessary to detect distortions of the solar- 
neutrino spectrum 

The P decay of *B populates a broad distribution of 
excitation energies in '^Be. Subsequently, the excited *Be 
nucleus breaks up into two a particles. The ^B neutrino 
spectrum can be deduced from the a-a relative-energy 
distribution, which can be determined in a laboratory 
experiment. Our collaboration has used two different 
experimental techniques to determine the a-a relative- 
energy distribution 0, Q . The two techniques give consis- 
tent distributions, which disagree slightly with the distri- 
butions determined in two previous experiments by Bhat- 
tacharya et al. [3| and Winter et al. [1, 111- In fact, the 
shapes of the former and the latest distributions are in 
reasonable agreement, but the maxima of the distribu- 
tions deviate by approximately 20 keV, which results in 
a 4% difference on the ^B neutrino spectrum at high en- 
ergies [El, > 14 MeV). This uncertainty complicates the 
attempts being made to detect neutrinos from the fusion 
of ^He+p in the Sun, and the attempts being made to 
detect the diffuse supernova neutrino background 0] ■ 

The purpose of the present Rapid Communication is 
to identify and quantify systematic effects not accounted 
for in the two former measurements of the a-a relative- 
energy distribution, which can explain the 20 keV de- 
viation relative to the two latest measurements per- 
formed by our collaboration, thus settling the current 
dispute 

Response of Si detectors to different ions. Both the 
former and the latest measurements were made with Si 
detectors, and in all four cases the energy calibration of 
these detectors was the main source of uncertainty. Typ- 



ically, the calibration of a Si detector consists in estab- 
lishing a linear relationship between the observed pulse 
height, X, and the deposited energy, E. However, the 
physical quantity that is directly related to the observed 
pulse height is not the deposited energy, but the number 
of collected electron-hole pairs, N. Therefore, it is more 
correct to assume a linear relationship between x and N, 
i.e.^ 

x = aN + b. (1) 

The energy deposited in the detector by a charged parti- 
cle can be divided into two parts: the part that goes into 
electron- hole pair formation (ionizing), which we shall 
denote by Ei, and the part that does not (non- ionizing), 
which wc shall denote by i.e. E = Ei + E^. Denoting 
the electron-hole pair collection efficiency by 77 and the 
energy required to create an electron-hole pair by e, we 
have 

N^^-^^liE-E^^). 

Introducing k = 7]/e, this simplifies to 

N^k{E~ En) . (2) 

Finally, combining Eqs. ([T|) and we find 

x = aK{E- E,,) + b. (3) 

Note that En depends on the charge and the mass of the 
particle and is a function of the deposited energy, E. Ide- 
ally, the parameter k would be independent of the charge 
and the mass and would be constant with E. However, 
the experimental study of Lennard et al. [Toj shows that 
K has a weak dependence on these quantities. (The in- 
terpretation of this effect suggested by Lennard et al. is 
that e decreases with increasing stopping power, but this 
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^ To simplify the discussion, we neglect non-linearities in the elec- 
tronic processing of the signal (amplification, shaping, and digital 
conversion). 
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is irrelevant to the present analysis.) Finally, the cali- 
bration coefEcients a and h are strictly independent of 
charge, mass, and E] they are determined by the exper- 
imental apparatus, i.e. the detector and the associated 
electronics. 

The experiment of Winter et al. In the experiment of 
Winter et al. 0, 0, the ions were implanted in a Si 
detector and the linearly added pulse heights of the two 
a particles were measured. The a particles have nearly 
equal energies, i.e. Ea^ = E^^ = ^E. The number of 
collected electron-hole pairs created by each a particle is 
given by Eq. The total number of collected electron- 
hole pairs is simply twice this number, i.e. 



Nsji^K^[E-2EnAhE)]' 



(4) 



The energy calibration was performed by implanting the 
/3-delayed a emitter ^°Na in the detector. The linearly 
added pulse heights of the a particle and the ^^0 ion 
were measured. In this case, the energy is shared approx- 



imately in ratios of 4:1, i.e. E^ = -^E and -Bieo 



and the total number of collected electron-hole pairs is 
given by 



N'. 



20 Na 



(5) 



In general, iVsg ^ N2o^g_, i.e. for the same deposited 
energy, E, different numbers of electron-hole pairs are 
collected and hence different pulse heights are measured. 
This is the basic reason why a calibration based on ^°Na 
data cannot be applied to ®B data without proper correc- 
tion. Note that Eqs. (|4|) and ([5|) do not include the small, 
but significant, contribution of the j3 particles, which we 
are not concerned with here. It follows from Eqs. (|4]) 
and ^ that the necessary correction for the different 
response of Si detectors to a particles and ^^O ions, is 
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FIG. 1: Reanalysis of the data of Lennard et al. [lol ]. Ewq is 
the energy deposited in the detector and 5Ewq is defined in 
Eq. ([7|. The dashed line shows the value adopted by us for 
5Ei6Q in the energy range 0.5-1.2 MeV. 



(pulse height) were determined experimentally at 13 en- 
ergies between 0.53 and 2.96 MeV. A linear fit to the 13 
data points (Ea — E^.a, ^a) gives^ 



0.1843(2) keV" 



0.2(3), 



with x^/dof — 1.3. Next, we consider the ^^O data. Cor- 
responding values of Ewq and xwq were determined ex- 
perimentally at four energies between 0.73 and 2.52 MeV. 
Using these data and the values of an^ and b determined 
above, we evaluate the following quantity. 



SEieci ~ El6r\ — 



X16Q 



aKo 



(7) 



AE = {Ns-s - 7V20Na 



IE 



K16Q 



He ~ E,,.eo{lE)] 



En,a{TE) — 2En^ai^E) 



(6) 



We will now determine the magnitude of AE using the 
experimental data of Lennard et al. [13] ■ In this experi- 
ment, the response of a Si detector to ^H, ^He, ^Li and 
^^O ions was studied. The energy of the incident ion {Ek 
in the notation of Lennard et al.) was determined by the 
time-of-flight technique. The energy loss in the detec- 
tor window {AEvj in the notation of Lennard et al.) was 
determined from the rate of change of pulse height with 
angle. The energy deposited in the detector could then be 
determined as E = Ek — AE^ . Finally, the non- ionizing 
energy deposition (i?n in the present notation and AEn 
in the notation of Lennard et al.) was estimated from 
theory. First, we use the a-particle data of Lennard et 
al. to calibrate the Si detector used in their experiment, 
i.e. we determine the parameters and b in Eq. ©. 
Corresponding values of E^ (deposited energy) and Xa 



which is shown in Fig. [TJ We conclude that SEwq is 
approximately constant in the energy range of interest 
with regard to the ^°Na calibration, Eibq = 0.5-1.2 MeV, 
and in this range SEwq = 66(4) keV. Substituting Eq. ([3]) 
for xwQ in Eq. ([7]), we obtain 



SE 



160 



K16Q 

El6Q — [Ei6q 



E, 



(Ev 



which we recognize as the first two terms in Eq. ([6]) with 
■^E = Ei6Q. Interpolating the theoretical calculations of 
Lennard et al., we determine the two last terms of Eq. ([6]) 
to be En,a{^E) = 9.1 keV and E^^^ci^E) = 9.8 keV 
with negligible uncertainties [l3|, for a total energy of 
E = 3.0 MeV corresponding to the maximum of the a-a 
relative-energy distribution in the /3 decay of ^B. We thus 



^ The data is consistent with an energy-independent Ka. In con- 
trast, we make no assumptions about kibq- 
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TABLE I: Maxima of the Be excitation-energy distribution 
obtained in the four most recent studies of the /3 decay of *B. 
The ^Be excitation energy, Ex, and the a-a relative energy, 
Ea-ct, are related as Ex = Ea-a — Egs, where Egs ~ 92 keV 
is the energy of the ground state of *Be relative to the 2a 
threshold. The deviation is computed with respect to the 
average of Refs. ^ and 'a']. 



Study 


Maximum 
(keV) 


la uncertainty 
(keV) 


Deviation 
(keV) 


Winter et al. \§] 


2 943 


9 


21 


Bhattacharya et al. [4] 


2 939 


5 


17 


Kirsebom et al. [2] 


2 921 


5 


-1 


Roger et al. [3] 


2 923" 


6 


1 



"This is the updated value given in Ref. 01 . 



conclude that 

AE = SEieo + ^„,«(|£;) - 2£;„,,(i£;) 

= 66(4) keV + 9.8 kcV - 2 x 9.1 keV 

= 58(4)keV. (8) 

This correction is approximately 13 kcV larger than the 
correction applied by Winter et al, which is said to be 
"40-50 keV for the various alpha lines, with an uncer- 
tainty of 5 keV" d il, citing Lennard et al. The 
sign and the magnitude of this difference is such that it 
can account for the 20 keV shift of the distribution of 
Winter et al. relative to our distributions, cf. Table [J 
We thus conclude that a misinterpretation of the data of 
Lennard et al. on the part of Winter et al. is the likely 
reason for the discrepancy between their result and ours. 

^°Na was used for energy calibration in both of our 
experiments, but only the implantation technique 0] is 
sensitive to the different response of Si detectors to a par- 
ticles and ions. The internal consistency of our ex- 
periments may thus been seen as a verification of Eq. ([8|). 

Winter et al. also used an external ^^^Th a source to 
provide calibration points at higher energies (5-9 MeV). 
Corrections were made for the energy loss of the a parti- 
cles in the source, but it is unclear whether the different 
implantation depths of the a-unstable daughter nuclei in 
the ^^^Th decay series were taken into account. As dis- 
cussed in Section III of Ref. , neglect of this effect will 
distort the energy calibration. 

Reanalysis of ^^B, ^^N, "^^Na decay data. To further 
support our case, wc have reanalyzed data from an ex- 
periment conducted by our collaboration where ^^B. -'^^N, 
and 2"Na ions were implanted in a finely segmented Si de- 
tector and the linearly added pulse heights of the charged 
decay products were measured [ll|. The same detector 
was later used in one of our experiments Q. The ^^B, 
^^N, and ^°Na energy spectra are shown in Figs. [2] (a) 
and (b). The two most prominent peaks in the ^''Na spcc- 




FIG. 2: (Color online) Decay energy-spectra of ^"Na ions (a) 
and ^^B and ^^N ions (b) implanted in a finely segmented Si 
detector [ij]. The peaks are labeled by the excitation energy 
in the daughter nucleus in MeV. 



trum have been used to calibrate the energy axis, taking 
into account the additional energy deposited by the _B 
particles, as determined with a GEANT4 simulation 0], 
and taking into account the /3-decay recoil energy, which 
is relatively small. The (3 decays of -'^^B and "'^^N pop- 
ulate excited states in ^^C, which breakup into three a 
particles. For the purpose of estimating the non-ionizing 
energy loss, £'n,Q, we assume that each a particle gets 
one-third of the breakup energy. This is a reasonable as- 
sumption for the two states in ^^C under consideration. 
Since E'n^Q only varies slowly with energy, about 0.9 keV 
per MeV, the error introduced by this approximation is, 
in any case, very small. In analogy to Eq. ([5]), we can 
then write down an expression for the error that results 
if the ^°Na calibration is applied to the ^^B and ^^N data 
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without correction for the different response of Si detec- 
tors to a particles and ^^O ions, 

= 66(4) keV + E^^^{^E) - iE^,c.{\E) . (9) 

The data analysis is summarized in Table [iTl As seen in 
Fig. [2] (b), the /3-decays of ^^B and ^^N populate two 
narrow states in ^^C, the 0+ state at 7.654 MeV and the 
1+ state at 12.710 MeV. When the additional energy de- 
posited by the /? particles and the /3-decay recoil energy 
are taken into account, we find that the peaks are located 
above their nominal energies by 54(4) keV and 49(9) keV 
for and by 56(4) keV and 61(7) keV for ^^N. These 
numbers are in good agreement with the deviations ex- 
pected based on Eq. 53(4) keV and 48(4) keV, thus 
lending further support to our case. 

The experiment of Bhattacharya et al. In the exper- 
iment of Bhattacharya et al. Q, the ^B ions were im- 
planted in a thin carbon foil viewed by an arrangement 
of Si detectors. Before the a particles can deposit their 
energy in the sensitive volume of the detectors, they must 
exit the foil and traverse the dead layer of the detec- 
tor, which leads to a small, but significant, energy loss. 
Interactions with electrons are frequent and lead to an 
approximately Gaussian energy-loss distribution, assum- 
ing that the layers traversed are not too thin, nor too 
thick [3: mi- Interactions with nuclei (Coulomb scat- 
tering) are far less frequent and do not contribute to the 
mean energy loss. However, because individual scatter- 
ing events in this case can result in a large energy loss, 
the energy-loss distribution acquires a pronounced high- 
energy tail, which is seen as low-energy tail in the re- 
sponse function of the detector. The magnitude of the 
tail is very sensitive to the thickness traversed [l^. In 
one of the two ^B experiments performed by our collab- 
oration 0, it was indeed found that the response func- 
tion determined with a ^^^Am a source had a more pro- 
nounced low-energy tail than the response function de- 
termined with ^'^Na implanted in a thin carbon foil. 

The effect of the low-energy tail is to shift the appar- 
ent position of the maximum of the a-particle energy 
spectrum toward lower energies. Bhattacharya et al. |j| 



used a ^ "Gd a source to determine the response function 
of their detectors, which could have led them to overes- 
timate the magnitude of the low-energy tail. This, in 
turn, would have led them to conclude that the actual 
position of the maximum of the a-particle energy spec- 
trum is higher in energy than it really is. If the difference 
were the same as in our experiment 0, the error would 
be about 8 keV, leading to a 2 x 8 = 16 keV error on 
the a-a relative energy. The sign and the magnitude of 
this error is such that it can account for the 20 keV shift 
relative to our distributions. In our experiment, the re- 
sponse function was extracted from the ^^Na data. Since 
the 2°Na ions and the ®B ions were implanted in the same 
carbon foil with similar depth profiles, the low-energy tail 
should be very similar in the two cases. 

Finally, we note that Bhattacharya et al. determine 
the peak of the a-a relative-energy distribution to be 
slightly wider than we do and Winter et al. do, suggesting 
that experimental broadening effects were not completely 
accounted for in their experiment. 

Summary and conclusion. Our collaboration has de- 
termined the a-a relative-energy distribution in the /3 
decay of ®B using two different experimental techniques, 
which yield consistent distributions 0, 0] that, however, 
deviate by a small, but significant, amount with respect 
to two previous determinations by Bhattacharya et al. Q 
and Winter et al. [11,0]. In the present Rapid Communi- 
cation, we have identified possible sources of error in the 
energy calibrations of the previous experiments that may 
account for the discrepancy. We suggest that Winter et 
al. did not correctly account for the different response 
of Si detectors to a particles and ^^O ions. We have 
presented empirical evidence supporting our suggestion 
and estimated the associated error to be approximately 
-1-13 keV. We suspect that Bhattacharya et al. may have 
overestimated the magnitude of the low-energy tail of 
their response function and we have argued that the as- 
sociated error could be as large as -1-16 keV. As an im- 
portant corollary, the ^B neutrino spectra deduced from 
the four experiments now agree within about 1% at high 
neutrino energies, whereas before they differed by 4%. 
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TABLE II: Reanalysis of the experimental data of Ref. [Tl|]. is the excitation energy of the daughter nucleus and E is 
the breakup energy, calculated as E = E^ - Sa for ^°Na and as E = E:^ — Ssa for ^^B and ^^N, where 5"^ = 4 729.8 keV 
is the a + ^^O threshold energy in ^^Ne and 5*30, = 7 274.7 keV is the 3a threshold energy in ^^C; both are known to high 
precision. Ejs is the shift due to the energy deposited by the P particles, determined as the mean of the energy-loss distribution 
obtained with a GEANT4 simulation -Br is the /3-decay recoil energy averaged over lepton angles, cf. Section IV in Ref. Q]- 
E' — E + Efs + i?R is the breakup energy corrected for the P- and recoil shifts. -Bcxp is the mean energy of the peaks in 
Figs. [2] (a) and (b). The ^"Na peaks were used to calibrate the energy axis. ASoxp = i?cxp — E' is the deviation with respect 
to the corrected breakup energy. Finally, AE is the expected deviation based on Eq. (|9}. All energies are in keV. 



E:, [12, 


E 






E' 






AE 


20j^^ 7421.9(12) 
10 273.2(19) 


2 692.1(12) 
5 543.4(19) 


26(2) 
31(2) 


0.5 
0.2 


2 719(2) 
5 575(3) 








i2g/i2j^ 7654.2(15) 
12 710(6) 


379.5(15) 

5 435(6) 


28(2)/26(2) 
59(2)/30(2) 


1.2/2.7 
0.0/0.5 


409(3)/408(3) 
5 494(6)/5 466(6) 


462(3)/464(3) 
5 543(6)/5 527(3) 


54(4)/56(4) 
49(9)/61(7) 


53(4)" 
48(4) 



"An additional correction of +1.3 keV was necessary to account for the energy dependence of -Bn,a outside the energy range covered 
by the ^''Na calibration data. 
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